Evolution of photoemission spectral functions in doped transition metal oxides 
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We discuss the experimental photoemission and inverse photoemission of early transition metal 
oxides, in the light of the dynamical mean field theory of correlated electrons which becomes exact in 
the limit of infinite dimensions. We argue that a comprehensive description of the experimental data 
requires spatial inhomogeneities and present a calculation of the evolution of the spectral function in 
an inhomogenous system with various degrees of inhomogeneity. We also point out that comparaison 
of experimental results and large d calculations require that the degree of correlation and disorder 
is larger in the surface than in the bulk. 
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There has been a resurgence of interest in doped tran- 
sition metal oxides following the discovery of high Tc 
oxides ||]. Early transition metal oxides, containing Ti 
and V, have attracted particular attention in very recent 
times La^rSri-j^TiOa ||] is such an example where 

the carrier concentration can be varied in a controlled 
fashion. Varying x tunes the electron configuration con- 
tinuously between and 1, with LaTiOa {x — 0) being 
a Mott-Hubbard insulator with a 3d} configuration and 
SrTiOs {x — 1) being a band insulator with a 3dP configu- 
ration. For large values of x, the transport and magnetic 
properties as well as the experimentally observed spec- 
tral weights behave as a correlated doped Mott insulator. 
For small values of x 0), the system represents small 
amount of electron doping in a 3cP configuration, and 
thus, it should behave as a doped band insulator. Indeed 
transport and magnetic properties of such compounds 
exhibit essentially free-electron like behavior. 

On the thoeretical side, it has been shown that the 
Hubbard model in the limit of large dimensions offers a 
semiquantitative description of the evolution of the trans- 
port and thermodynamic data as a function of carrier 
concentration, in these systems. The consistent mod- 
elling of the photoemission spectroscopy, however, is more 
problematic and is the subject of this paper. 

We first summarize the main experimental facts that 
a theory should describe. 1) In a pioneering paper Fuji- 
mori et al. 1^ investigated the spectral function of sys- 
tems with a d^ configuration. They observed that there 
are two features in the UPS spectra, a low energy fea- 
ture near the Fermi edge, and an intermediate energy 
feature around 1.4 eV which we refer to as the coherent 
and the incoherent parts of the spectral function respec- 
tively. The metal-to-insulator transition at filling n = 1 
is driven by the transfer of spectral weight from the low 
energy feature to the incoherent part of the spectra. In- 
oue and collaborators confirmed this picture in a more 
careful study of the CaSrVOs system This behavior 
can be understood in terms of the solution of the infinite 
dimensional Hubbard model at half filling but the 
parameters required to fit the data produce specific heat 
and susceptibility enhancements which are inconsistent 
with the thermodynamic data in these systems. [ p| 

2) In a subsequent publication Fujimori and collabora- 
tors addressed [|| the doping dependence of the coherent 
and incoherent features in Laa;Sri_a;Ti03 with various 
values of x . They found that as one dopes the Mott in- 
sulator, the weight in the incoherent and in the coherent 
features grow. Furthermore, they observed that the inco- 
herent feature persists up to fairly large values of doping, 
when one would have expected nearly free electron behav- 
ior. Similar behavior was reported in the Yi_xCaxTiO^ 
system. The behavior at large doping (small electron 
number) has been confirmed in careful studies of the 
SrTi03-5 system. [| 

In contrast to the experiments described in the previ- 



ous paragrpah, no theoretical model ever came close to 
explaining these observations. In particular the large d 
Hubbard model which successfully accounted the trans- 
port and thermodynamic data in the doped Mott insu- 
lator, completely fails in reproducing the experimentally 
observed evolution of the spectral function as one varies 
the carrier concentration. 

3) Sarma and collaborators revisited the limit of small 
density of electrons. They performed photoemission 
and inverse photoemission spectroscopy on SrTiOs^^ 
where the oxygen deficiency 5 introduces electrons into 
the stoichiometric d° system. The photoemission part of 
the spectra confirmed the findings of Fujimori et al. with 
an incoherent feature around 1.4 eV and a coherent fea- 
ture at low frequencies, while the inverse photoemission 
spectra were in complete agreement with band structure 
calculations. 

We stress again that these results are very surprising, 
since in a dilute system the interactions are always renor- 
malized to small values so both photoemission and in- 
verse photoemission should be predicted correctly by the 
band structure calculations. 

4) The discrepancy between the parameters that are 
necessary to fit the bulk thermodynamical and transport 
properties on one hand and the highly surface sensitive 
photoemission spectroscopy on the other hand, which 
was alluded to earlier, suggests that the behavior near 
the surface may be different than in the bulk. While 
this issue is technically difficult to be probed on the di- 
lute d^ systems such as La^^Sri-a, TiOa or SrTiOg-a, it 
can be effectively addressed in the denser systems such 
as the LaajCai-xVOa system with electron configuration 
varying between 1 and 2. Sarma and coworkers Q per- 
formed a systematic study of the evolution of the spec- 
tral function as a function of photon energy and observed 
that the incoherent features get weaker while the coher- 
ent features become stronger as the probing depth of the 
spectroscopic technique increases by varying the photon 
energy suitably. 

Taking into account all these disparate experimental 
observations, we believe that to obtain a coherent pic- 
ture of the experimental situation in these systems one 
has to recognize, that 1) a minimal model of the spec- 
troscopic data in these systems has to recognize that the 
surface experiences a different degree of correlation and 
disorder than the bulk. 2) The ultra-violet photoemis- 
sion spectroscopy can only be modelled successfully by 
models that include spatial inhomogeneities. We argue 
that the existent data are good enough to rule out ho- 
mogeneous models but cannot determine the length scale 
over which the inhomogeneity takes place. 

In the framework of the infinite dimensional Hubbard 
model we present calculations of the evolution of the 
spectral function with increasing doping in an interact- 
ing disordered system and argue that they are in good 
agreement with all the experimental observations cited 
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above. Whereas substituted compounds have been ex- 
tensively investigated in recent times, the role of disorder 
has almost always been neglected; thus our findings have 
important implications for a large number of materials. 
The present work also brings into focus the different de- 
grees of disorder and correlations at the surface and in the 
bulk; while this is indeed to be expected and should be 
taken into account when interpreting experimental data, 
but has not received any attention before. 

In order to motivate our model we note that the Hub- 
bard model represents a homogeneous system, whereas in 
real systems electron (or hole) doping is brought about 
by chemical substitution; this invariably introduces dis- 
order or inhomogeneities at some length scale. There are 
several reasons why the strengths of the correlations and 
the disorder are different in the bulk than at the surface. 
Atoms at the surface experience a different coordination 
number than in the bulk, therefore the ratio of the inter- 
action strength to the average hopping matrix element is 
larger at the surface than in the bulk. The strength of 
the disorder is related to the spatial fluctuations in the 
effective potential experienced by an atom. In the limit 
of infinite coordination off-diagonal disorder becomes a 
non fluctuating quantity , i. e. it is self averaging. So the 
effective strength of off-diagonal disorder increases with 
decreasing effective dimensionality. 

In this paper we propose that it is essential to use a spa- 
tially inhomogeneous picture to model the surface of these 
systems. The basic idea is that in these systems the elec- 
trons experience different local environments. The an- 
gle integrated photoemission spectrum, as a local probe, 
captures a suitable average of the Greens function over 
the local environments. In principle the degree of inho- 
mogeneity can vary between two limits. If the system 
is disordered on a microscopic scale (termed microscopic 
inhomogeneity here), but homogeneous on a macroscopic 
scale, i.e. the impurities do not form clusters, a model 
of uncorrelated disorder is appropriate. If there is a ten- 
dency of the defects to form clusters, this may reinforce 
the natural tendency of strongly correlated systems to 
phase separate as pointed out in ref and a model of 
macroscopic inhomogeneities is more suitable. In this pa- 
per we use the dynamical mean field theory that becomes 
exact in the limit of infinite spatial coordination [p^ to 
make these ideas more quantitative. We study the evolu- 
tion of photoemission spectra in spatially inhomogenous 
systems which we model using diagonal disorder. Off 
diagonal disorder which would correspond more closely 
to the physical situation disappears in the limit of large 
dimensions, unless one introduces very special distribu- 
tions of hopping matrix elements . We find that the 
photoemission spectrum as a function of x to be very 
different from the predictions of the Hubbard model and 
very similar to what is observed experimentally. 

We study a Hubbard model with diagonal disorder: 



(ij) i « 

To get a bounded density of states and a well defined 
limit as the lattice coordination gets large, we take the 
hopping matrix elements tij to be defined on a Bethe 
lattice and scale them as The full bandwidth 2D — 
4t, is estimated to be 2.5 eV to fit our band structure 
calculations. The are random variables which take the 
values eA with probability Pa and es with probability 
Pb ; Pa and Pb are given by the composition, as 2S and 
(1 — 26) in SrTi03_5 and x and (1 — x) in LaajSri-^jTiOa 
respectively. Supercell calculations for SrTi02.875 and 
SrTi02.75 have shown that several bands are pulled down 
compared to SrTiOs due to the presence of vacancies 
by about 1.5 - 3.0 eV which provides an estimate for 
(es — e^)- We use (es — e^) ~ 2.4 eV; however, the 
qualitative features of the calculations do not depend on 
this exact value. Ua = Ub = 4.6 eV was estimated from 
earlier spectroscopic analysis. The chemical potential fi 
is chosen so that the number of d electrons is equal to 
X in La^jSr^^ — ^ 

TiOa or 26 in SrTiOs-^. The diagonal 
disorder introduces randomness that survives the limit of 
large dimensionality, and captures the basic physics that 
different sites experience different local enviroments. 

The mean field theory which is exact in the d ~^ oo 
limit determines the local Greens functions at the A and 
B sites as the local / electron Greens function of an An- 
derson impurity model 

k 

+ E ^'^("t/- + f^f^o) + eff+f. (1) 
fc 

Here the / level position e/ is given hy e a — fJ- and sb — fJ-, 
respectively. The hybridization function A(ia;„) 

A(^c.„) = E7— ^ (2) 

has to be determined self consistently from the local 
Greens functions via the equation: 

A(iLj„) = t^ Gav[ii^n) = {Pa GA{iuJn) + PB GB{iuJn)) 

It is to be noticed that the A and B sites are connected 
by the same hybridization function which depends on 
the local physics at both sites. Thus we have a coupled 
system of equations. To solve the mean field equations 
we use the recent extension of the IPT method to 
extract real frequency information. 

Figure (a) shows the total as well as the partial spec- 
tral functions for small filling (n = 0.2). Since >> e^, 
the particle concentration on the B sites is very small. 
Thus the corresponding spectral function is only slightly 
renormalized by interactions. It basically consists of a 
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semicircle around cb — fJ- with a small correlation induced 
peak at — /i + C^- On the A sites, in contrast, we find 
almost one particle per lattice site. Therefore ImGA(i^) 
exhibits features similar to a slightly doped Mott insula- 
tor, i. e. there are two Hubbard bands around eA — fJ- and 
eA — ^J■ + U plus a resonance peak at zero frequency. No- 
tice, that the lower Hubbard band of Ga(w) produces a 
correlation induced feature in the averaged Greens func- 
tion Gav{i^), as observed in the experiments. This max- 
imum grows with increasing doping as can be seen from 
the evolution of the spectral function shown in the left 
column of figure . In the limit n 1, there are only A 
sites left and the incoherent feature observed at negative 
frequencies becomes the lower Hubbard band of a pure 
Mott insulator. 

We now turn to the case where the inhomogeneities oc- 
cur on a macroscopic scale. In this scenario the density of 
electrons is inhomogeneous on a more macroscopic scale. 
A fraction /i is in an insulating phase with the local con- 
centration of the Mott insulating phase riins = 1 and a 
fraction 1 — /i with the metallic concentration n„iet- The 
total density n = finins-'r{y — fijnmet, and the measured 
Greens function is given by 

G = /iG„,,„ + (1 - /i)G„„,, (3) 

G„.„^ and G„^^j can again be calculated using the mean 
field theory now in the respective homogeneous phases. 
To explore this scenario we took Timet = -1 (which corre- 
sponds to a chemical potential /.t = —.78) and varied the 
density n by varying the volume fraction /i in equation 
. The total and the partial spectra for this case are 
shown in figure (b). Except for the absence of a res- 
onance peak at zero frequency the picture is similar to 
the case of microscopic disorder discussed above. In par- 
ticular, the lower Hubbard band of the insulating phase 
produces an incoherent feature at negative frequencies, 
which grows with increasing filling (see right column of 
figure ). The experimentally measured quantitiy is a con- 
volution of spectral distribution at negative (photoemis- 
sion) or positive (Bremsstrahlung isochromat) frequen- 
cies with a resolution function, (a) and (b) shows this 
convolution that can now be compared directly with the 
experimental data of Sarma et al. ||]. The data are in 
excellent agreement with the experimental spectra. For 
example, the calculation correctly predicts the existence 
of an intense incoherent feature at about 1.3 eV arising 
from the local Greens function for site A with a band-like 
coherent feature near Ep due to the local Greens function 
for site B. This is in perfect agreement with the experi- 
mentally observed features in the photoemission spectra 
of SrTi03-5 11 and La^Sri.^TiOa §. Moreover, it is 
clear from the calculated results that unoccupied parts of 
the spectral function will be entirely dominated by coher- 
ent band-like features arising from the site B as observed 
in the experiment ||^. 



On the other hand, we would like to stress that the 
experimentally measured specific heat, a bulk property, 
cannot be described using the parameters obtained from 
the (surface sensitive) photoemission experiments. Ex- 
perimental measurements Q show that 7 = C/T is en- 
hanced as the Mott insulator is approached with increas- 
ing filling. It is clear that this observation cannot be 
accounted for by phase separation, because the metallic 
phase is assumed to be a dilute system hardly showing 
any correlation eff'ects and the insulator does not con- 
tribute to the linear coefficient of the specific heat. ( we 
assume that the insulating state gives a negligible contri- 
bution to the specific heat ) The scenario of microscopic 
disorder in contrast vastly overestimates this quantity. 
This can be understood by the fact that the electrons 
get trapped on the A sites so that the local particle den- 
sity ua is close to one and much larger than the average 
particle density n (see above). As has been established 
in reference [0, the specific heat can, however, be de- 
scribed by the homogeneous Hubbard model. 

To summarize, in this paper we discuss a model for the 
high-energy spectroscopic data of compounds with close 
to 3(i" configuration. We presented explicit calculations 
of the spectral function of correlated electrons at different 
dopings in an inhomogeneous picture. We considered two 
cases; in one the inhomogeneity is present only on a mi- 
croscopic scale while in the other the system was assumed 
to be macroscopically inhomogeneous. Both models give 
similar angle integrated spectra which therefore cannot 
determine the length scale over which the system is in- 
homogeneous; however, the calculations are in excellent 
agreement with known experimental results for such and 
similar systems. The spectra of the inhomogeneous mod- 
els are qualitatively different from that of a homogeneous 
system, suggesting the need to go beyond the standard 
theories to describe spectral functions in a class of doped 
transition metal compounds. 
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FIGURES: 

• Fig. 1 Spectral functions of the disordered Hub- 
bard model for density n=.2, U=3.7 in the case of 
(a) microscopically and (b) macroscopically inho- 
mogeneous disorder. The energy scale is given in 
units of D (= 1.25 eV). 

• Fig. 2 Evolution of the spectral functions of the dis- 
ordered Hubbard model for U = 3.7, et — ea = 1.9, 
and increasing filling n; left column: homogeneous 
disorder, right column: inhomogeneous disorder 

• Fig. 3 Theoretical spectra of the disordered Hub- 
bard model for density n = .2, U = 3.7, and 
£(, — Eo = 1.9 in the case of (a) microscopic and (b) 
macroscopic inhomogeneities. The dashed and dot- 
ted curves indicate the partial contributions from 
the A and B sites. The occupied parts of the spec- 
tral functions are expanded for clarity, as shown in 
the figures. 
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